Anisotropic Surface Growth Model in Disordered Media 
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We introduce a self-organized surface growth model in 2+1 dimensions with anisotropic avalanche 
process, which is expected to be in the universality class of the anisotropic quenched Kardar-Parisi- 
Zhang equation with alternative signs of the nonlinear KPZ terms. It turns out that the surface 
height correlation functions in each direction scales distinctively. The anisotropic behavior is at- 
tributed to the asymmetric behavior of the quenched KPZ equation in 1+1 dimensions with respect 
to the sign of the nonlinear KPZ term. 



-a 
c 
o 

> 

in 
o 

o 

\o 

On 

-i— > 
c3 



i 

C 

O 

o 



X 



PACS numbers: 68.35.Fx, 05.40. +j, 64.60.Ht 



Subject of the pinning-depinning (PD) transition by 
external driving force has been of much interest recently. 
The problems of the interface growth in porous media 
under external pressure jl ] , the dynamics of domain wall 
under random field the dynamics of charge density 
wave under external field ||, and the vortex motion in 
superconductors under external current [QJ|] are typical 
examples. In the PD transition, there exists a critical 
value F c of the driving force F, such that when F < F c , 
interface (or charge, or vortex) is pinned by disorder, 
while for F > F c , it moves with constant velocity v. The 
velocity v plays the role of order parameter in the PD 
transition which behaves as 



(F - F c ) 6 



(1) 



Recently, several stochastic models for the PD tran- 
sition of interface growth in disordered media have been 
introduced It is believed that the models in 1+1 di- 
mensions, are described by the quenched Kardar-Parisi- 
Zhang (QKPZ) equation 

d t h = vdlh+hd x hf + F + j 1 {x,h), (2) 



where the noise r\ depends on position x and height h with 
the properties of (rj(x,h)) — and (r)(x, h)rj(x', h')) = 
2D5(x - x')S(h - ti). The QKPZ equation exhibits 
the PD transition at F c . The surface at F c can be de- 
scribed by the directed percolation (DP) cluster spanned 
perpendicularly to the surface growth direction in 1+1 
dimensions. The roughness exponent a of the interface is 
given as the ratio of the correlation length exponents of 
the DP cluster in perpendicular and parallel directions, 
that is a = v±/v\\ ~ 0.63. 

The origin of the nonlinear term in the QKPZ equation 
is different from that of the thermal KPZ equation with 
the noise r/(x, t) For the quenched case, the nonlinear 
term is induced by the anisotropic nature of disordered 
media, while for the thermal case, it is induced by lateral 



growth, and thus the coefficient A is proportional to the 
velocity of the interface, which vanishes at the thresh- 
old of the PD transition. For surfaces belonging to the 
DP universality class, the positive nonlinear KPZ term 

is induced under coarse graining of the quenched random 

r /2 l /2 . 

force with amplitudes A h and A x in /i-direction and 

in x-direction, respectively, when > A x . On the other 
hand, one may consider the case in 2+1 dimensions that 
the amplitudes of random force are anisotropic on sub- 
strate, that is, Ah > Am in one direction of the substrate 
and Ah < Aj_ in the other. In such a case, following [8], 
the coarse-grained Langevin equation is expected to take 
the form of the anisotropic QKPZ (AQKPZ) equation 
given by 



dth = v\\d\\ 2 h 



d L 2 h 



+F + V (r,h), (3) 



with Ay > and A i < 0. 



In this letter, we study the surface of Eq. (3) by intro- 
ducing a self-organized stochastic model. The universal- 
ity class of the stochastic model is checked by comparing 
surface properties with those obtained from direct nu- 
merical integration of Eq. (3). As shown in Fig. 1, the 
surface of Eq. (3) at F c forms the shape of a mountain 
range with steep inclination in one direction, whereas 
it is gently sloping in the other direction. Accordingly, 
the roughness exponents of the height-height correlation 
functions in each direction scale distinctively, which leads 
to a new universality class. This result is remarkable as 
compared with the case of the thermal noise. For the 
thermal case, the height-height correlation function is 
isotropic, and the anisotropic KPZ equation in 2+1 di- 
mensions renormalizes into the Edwards- Wilkinson equa- 
tion p. This is because for the case described by the 
KPZ equation with thermal noise, the sign of the non- 
linear terms is irrelevant in determining the universality 
class p0| , and the two nonlinear terms with different 
signs are canceled out effectively [[ll). However, for the 



quenched case, the anisotropic surface morphologies of 
Fig. 1 imply that such cancellation does not occur and 
the QKPZ equation is asymmetric with respect to the 
sign of the nonlinear KPZ terms. Accordingly, in this 
letter, we also study the QKPZ equation with negative 
A in 1+1 dimensions, which enables one to understand 
the anisotropic nature of the AQKPZ equation in 2+1 
dimensions. 

First, we study the QKPZ equation with A < in 1+1 
dimensions by direct numerical integration with the dis- 
cretized version, 

h(x, t+At) = h(x, t) + At{h(x - 1, t) + h(x + 1, t) 

l,t)) 2 +F} 



-2h{x, t) + -(h(x + 1, t) - h{x 
8 

+(At) 2 / 3 e(x, [h(x,t)]), 



(4) 



where [• • •] denotes the integer part, and £ is uniformly 
distributed in |]. The prefactor (At) 2/3 of the 

noise term arises from approximately coarse-graining the 
noise r)(x,h) during the time interval At. Numerical 
integration using Eq. (4) for A = 1 and At — 0.01 yields, 
even for a modest system size of L = 10 3 , the roughness 
exponent a ~ 0.63, which is consistent with the value 
of the DP universality. Note that the use of usual pref- 
actor (At) requires much larger computational cost to 
obtain the DP value of a ]l2|,|3|. Details for derivations 
of Eq. (4) and result of the direct numerical integration 
for A > will be published elsewhere fll4| . 

For the case of A < 0, we performed the numerical 
integration with A = — 1 and At = 0.01 for convenience. 
Fig. 2 shows typical surface configurations evolved tem- 
porally at the PD transition point, F c m 1.98, which 
exhibits the shape of a mountain with flat inclination in 
the pinned state. The surface in the pinned state looks 
similar to that of the model A by Sneppen Jl5| , and the 
shape of the surface determines the roughness exponent 
to be a = 1. Note that the RSOS restriction before 
deposition in the model A of Sneppen does not allow the 
particle to deposit in every site, which makes the growth 
velocity reduced, and results in A < jL(|. Accordingly, 
it is reasonable to have the surface morphology as shown 
in Fig. 2 for the QKPZ equation with negative A, which 
is different from the one for A > 0. In addition, we 
also examined the noise distribution on perimeter sites. 
It reveals that the pinning is caused by relatively large 
pinning strengths around the site where the height is 
minimum. The flatness on inclination makes the term 
(d x h) 2 large, which with the large pinning strengths 
compensates the external driving force. Thus the growth 
velocity becomes zero, and the surface is pinned. Since 
the surface pinning is caused mainly by the barrier of the 
pinning strengths around the site of minimum height, the 
growth velocity exhibits a sudden jump as the barrier is 
overcome by increasing the external driving force F. Ac- 
cordingly, the PD transition is of first-order as depicted 



in Fig. 3. On the other hand, for A > 0, pinned sites are 
not localized, but scattered, so that the PD transition is 
continuous. Therefore the surface of the QKPZ equation 
in 1+1 dimensions is asymmetric with respect to the sign 
of the nonlinear term. 



Next, we study the AQKPZ equation, Eq. (3), in 2+1 
dimensions by introducing a stochastic model. It is a nat- 
ural extension of the RSOS model introduced previously 
by the current authors to study the anisotropic thermal 
KPZ equation fTi| . The stochastic model is based on a 
combination of the Sneppen's model A for A < and 
the model B for A > which is realized by assigning an 
anisotropic avalanche process. The advantage of study- 
ing such stochastic model is twofold; first there is no need 
for fine-tuning of F to get the critical state and second, 
asymptotic states can be readily reached for small system 
size. The model is defined on the checkerboard lattice, 
square lattice rotated by 45°. Initially we begin with a 
flat surface characterized by the heights on one sublat- 
tice and 1 on the other (see Fig. 4). Random numbers 
are assigned to each site. At each time step, selected is 
the site with minimum random number among the sites 
of which two nearest neighbors at (i + h,j — h) 
and (£— h,j— \) are higher. The site is updated by in- 
creasing the height by 2. Next, the anisotropic avalanche 
process may occur on the neighboring sites, (i + i, j + |) 
and [i — h,j + \)- If their height is lower by 3 than 
that of then the height is increased by 2. The 

avalanche rule is then applied successively to next rows 
in j-direction until there is no change. The sites with 
increased height are updated by new random numbers. 
The avalanche direction, j-direction, corresponds to the 
r|| -direction in Eq. (3). The anisotropic avalanche pro- 
cess along positive j-direction is an interesting aspect 
of our model, which is a generalization of the model 
by Maslov and Zhang jl^] into 2+1 dimensions. Such 
anisotropic avalanche process is to be distinguished from 
the isotropic avalanche process on tilted substrates, in 
which the roughness exponent along the tilt direction 
(rx-direction) is 1/3 ||. Using the tilt argument, it can 
be shown that our model includes alternative signs of the 
nonlinear terms, that is, A|| > and A^ < 17|. A typi- 
cal surface morphology is shown in Fig. la. We measured 
the roughness exponents for the height-height correlation 



functions, C\\ (r\\) 
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and Cj_ (rj_ ) ~ r ±~ L ■ The roughness exponents for each 
direction are obtained as a\\ = 0.25(1) and a± = 0.75(1) 
as shown in Fig. 5. We also measured the height fluctua- 
tion width, W 2 = J2A h r ~ hf ~ L2a for t > L z and 
~ t 2f3 for t <^ L z , where the exponents a, f3, and z are 
the roughness, the growth, and the dynamic exponents, 
respectively, and L is system size. It is obtained that 
a = 0.87(1), j3 = 0.80(1), and z is given as z — a/(3. 
The values of the exponents are different from those of 
the isotropic case in 2+1 dimensions where a w 0.48 



2 



and (3 « 0.41 We also measured the avalanche size 
distribution, P(s) ~ s~ T . The exponent r is obtained as 
« 1.35(3). 



In order to check if the stochastic model reduces to 
the AQKPZ universality, we considered the surface of 
Eq. (3) by carrying out the direct numerical integration 
using the two dimensional version of Eq. (4). We used 
the numerical values of At = 0.01, A|| = 1, and A^ = —1 
for convenience. Fig. lb shows the surface morphol- 
ogy obtained by the direct numerical integration at the 
threshold of the PD transition, F c « 0.50, which looks 
similar to the one in Fig. la. However, we could not 
measure the roughness exponents ay and a± precisely, 
because their precise measurement requires relatively 
large system size and huge computing times. Neverthe- 
less, since the morphologies of Fig. la and lb are similar 
to each other and that of An > and Aj^ < can be 
proven using the tilt argument for the stochastic model 
jl7j , we believe that the stochastic model belongs to the 
AQKPZ universality. The PD transition turns out to be 
continuous as depicted in Fig. 6. The velocity exponent 
9 defined in Eq. (1) is obtained as 9 = 0.9(1) which is 
somewhat larger than 9 ~ 0.8 for the isotropic case |p^| . 
The numerical results, t sa 1.35, z « 1.09, and 9 « 0.9 
seem to represent the characteristics of the self-organized 
critical depinning transition p8| , p0| , but the relation of 
those exponents to the anisotropic roughness exponents 
in 2+1 dimensions is not clear yet. Further study is 
required about this point. We have also examined the 
surface in moving state, F > F c . In this regime, the sur- 
face is no longer anisotropic, and reduces to the AKPZ 
equation with thermal noise for F ^> F c . 



In summary, we have studied the AQKPZ equation 
with alternative signs of the nonlinear terms in 2+1 di- 
mensions, and found that it leads to a new universality 
class. The surface exhibits anisotropic scaling behavior, 
which is due to the asymmetric behavior of the QKPZ 
equation with respect to the sign of the nonlinear term. 
The QKPZ equation with A < in 1+1 dimensions has 
also been studied. We have obtained that the surface 
forms the shape of mountain with flat inclination and 
the PD transition is of first order. Since the anisotropic 
KPZ equation with thermal noise has been applied to 
the flux line dynamics |f2l| , the AQKPZ equation consid- 
ered in this letter may also be relevant to the flux line 
depinning problem in disordered media. Further details 
will be published elsewhere. 
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FIG. 1. Typical surface configuration of the AQKPZ equa- 
tion generated by (a) the stochastic model and (b) direct nu- 
merical integration at the transition point. 
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FIG. 4. The configuration of a substrate with one step 
along a row for linear size L — 4. The avalanche process 
occurs in j-direction. 
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FIG. 2. Temporally evolved surface configurations of the 
QKPZ equation with A < in 1+1 dimensions at the pin- 
ning-depinning transition point. Successive height profiles are 
shown at constant time intervals. 
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FIG. 5. The height-height correlation functions in paral- 
lel (lower data) and perpendicular (upper ones) directions for 
the AQKPZ stochastic model in 2+1 dimensions. 
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FIG. 3. The growth velocity versus force for the QKPZ 
equation with negative A in 1+1 dimensions. 
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FIG. 6. The growth velocity versus force for the AQKPZ 
equation in 2+1 dimensions. 
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